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We measure the lifetime r^o and the B^-B'^ oscillation frequency Am^ with a sample of 
approximately 14,000 exclusively reconstructed -6° — > D*^£'^ui signal events, selected from 23 
million BB pairs recorded at the T(4S') resonance with the BAR4R detector at the Stanford Linear 
Accelerator Center. The decay position of the other B is determined with the remaining tracks in 
the event, and its 6-quark flavor at the time of decay is determined with a tagging algorithm that 
exploits tlic correlation between the flavor of the b-quark and the charges of its decay products. 
The lifetime and oscillation frequency are measured simultaneously with an unbinned maximum- 
likelihood fit that uses, for each event, the measured difference in decay times of the two B mesons 
(Af), the calculated uncertainty on At, the signal and background probabilities, and 6-quark tagging 
information for the other B. The results are 

Tgo = (1.523io:o23 ± 0.022) ps 

and 

Amd = (0.492 ± 0.018 ± 0.013) ps"\ 
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The statistical correlation coefficient between Tga and Amd is —0.22. 
PACS numbers: 13.25.Hw, 12.15.Hh, 14.40.Nd, ll.30.Er 



I. INTRODUCTION AND ANALYSIS 
OVERVIEW 

The time evolution of mesons is governed by the 
overall decay rate l/Tgo and the B^-B^ oscillation fre- 
quency Am^. The phenomenon of particle-antiparticle 
oscillations or "mixing" has been observed in neutral 
mesons containing a down quark and either a strange 
quark [K mesons) or a bottom quark [B mesons) . 
In the Standard Model of particle physics, mixing is 
the result of second-order charged weak interactions in- 
volving box diagrams containing virtual quarks with 
charge 2/3. In B mixing, the diagrams containing the 
top quark dominate due to the large mass of the top 
quark. Therefore, the mixing frequency is sensitive to 
the Cabibbo-Kobayashi-Maskawa quark-mixing matrix 
element Vtd f^- In the neutral K meson system, mix- 
ing also has contributions from real intermediate states 
accessible to both a K'^ and a meson. Real interme- 
diate states lead to a difference in the decay rate for the 
two mass eigenstates of the neutral meson system. For 
the B system, the decay rate difference is expected to 
be of 0(10-2-10-3) times smaller [g than the average 
decay rate and the mixing frequency, and is ignored in 
this analysis. 

We present a simultaneous measurement of the B^ 
lifetime and oscillation frequency based on a sample of 
approximately 141)00 exclusively reconstructed B° — > 
D*~ t^vi decays selected from a sample of 23 million 
BB events recorded at the T(45') resonance with the 
BABAR detector Q at the Stanford Linear Accelerator 
Center, in 1999-2000. In this experiment, 9-GeV elec- 
trons and 3.1-GeV positrons, circulating in the PEP-II 
storage ring annihilate to produce BB pairs moving 
along the e'~ beam direction (z axis) with a Lorentz boost 
of (3^ = 0.55, allowing a measurement of the proper time 
difference between the two B decays. At. 

The decay-time difference At between two neutral B 
mesons produced in a coherent P-wave state in an T(4S') 
event is governed by the probabilities to observe an un- 
mixed event, 

P(B°B° ^ P°B°) cx e-l^*l/"««(l + cos AmdAt), (1) 
or a mixed event, 

P(B°B° ^ B"P" or P"P°)cxe-l'^*l/^i50(l-cosAm<jAi). 

(2) 

Therefore, if we measure At and identify the 6-quark fla- 
vor of both B mesons at their time of decay, we can 
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extract t^o and Atoc;. In this analysis, one B is recon- 
structed in the mode P° f£, which has a mea- 
sured branching fraction of (4.60 ± 0.21)% d). Although 
the neutrino cannot be detected, the requirement of a re- 
constructed D*^ — !■ D^iT~ decay and a high-momentum 
lepton satisfying kinematic constraints consistent with 
the decay P° — > Z3*-£+j/^ allows the isolation of a sig- 
nal sample with (65-89)% purity, depending on the D° 
decay mode and whether the lepton candidate is an elec- 
tron or a muon. The charges of the final-state parti- 
cles identify the meson as a P° or a B'^ . The remaining 
charged particles in the event, which originate from the 
other B (referred to as Ptag), are used to identify, or 
"tag", its flavor as a P'^ or a P°. The time difference 
At = to* I — ttag ~ l^z/ fi^c is determined from the sepa- 
ration Az of the decay vertices for the D*~£^ candidate 
and the tagging P along the boost direction. The average 
separation is about 250 /im. 

The oscillation frequency and the average lifetime of 
the neutral P meson are determined simultaneously with 
an unbinned maximum-likelihood fit to the measured At 
distributions of events that are classified as mixed and 
unmixed. This is in contrast to most published measure- 
ments ^ in which only Tgo is measured, or Anid is 
measured with Tgo fixed to the world average. There are 
several reasons to measure the lifetime and oscillation fre- 
quency simultaneously. The statistical precision of this 
measurement for both Tgo and ArUd is comparable to the 
uncertainty on the world average; hence, it is appropriate 
to measure both quantities rather than fixing the lifetime 
to the world average. Since mixed and unmixed events 
have different At distributions, the separation of mixed 
and unmixed events gives greater sensitivity to the At 
resolution function; as a result, the statistical uncertainty 
of Tgo is improved by approximately 15% |jl^. Also, since 
B^ B^ events do not mix, we can use the At distributions 
for mixed and unmixed events to help discriminate be- 
tween B'^B^ signal events and B'^B~ background events 
in the lifetime and mixing measurement. 

There are three main experimental complications that 
affect the At distributions given in Eqs. and |. First, 
the tagging algorithm, which classifies events into cate- 
gories c depending on the source of the available tagging 
information, incorrectly identifies the flavor of Ptag with 
a probability Wc with a consequent reduction of the ob- 
served amplitude for the mixing oscillation by a factor 
(1 — 2'Wc)- Second, the resolution for At is comparable to 
the lifetime and must be well understood. The probabil- 
ity density functions (PDF's) for the unmixed (-I-) and 
mixed (— ) signal events can be expressed as the convo- 
lution of the underlying Attruc distribution for tagging 



7 



category c, 



lAttrucl/l-gO 



- [1 ± (1 - 2wc) cos AmrfAttruc] , (3) 



with a resolution function that depends on a set of pa- 
rameters determined from the data. A final complication 
is that the sample of selected B'^ D*~£'^vg candidates 
includes several types of background for which the At 
distributions must be determined. 

To characterize the backgrounds, we select control 
samples of events enhanced in each type of background 
and determine the signal and the background probabili- 
ties for each event in the signal samples and the back- 
ground control samples as described in Sec. [V. The 



measurement of Az and the determination of At and 
the uncertainty on At ((lAt) for each event is discussed 
in S ec. M . The 6-qu ark tagging algorithm is described in 
Sec. |Vlfln Sec. VII, we describe the unbinned maximum- 
likelihood fit. The physics model and At resolution func- 
tion used to describe the measured At distribution for 
signal are given in Sec. VIII. A combination of Monte 



Carlo simulation and data samples are used to deter- 
mine the parameterization of the PDF's to describe the 
At distribution for each type of background, as described 
in Sec. IX. The likelihood is maximized in a simultaneous 



fit to the signal and background control samples to ex- 
tract the lifetime r^o, the mixing frequency Am^, the 
mistag probabilities Wc^ the signal At resolution param- 
eters qc^ the background At model parameters, and the 
fraction of 5+ D*^ l^v^X decays in the signal sample. 
The results of the fit are given in Sec. Cross-checks 
are described in Sec . |X] and systematic uncertainties are 
summarized in Sec. KIl . 



II. THE BABMl DETECTOR 

The BABAR detector is described in detail elsewhere [|| . 
The momenta of charged particles are measured with 
a combination of a five-layer silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH) in a 1.5- 
T solenoidal magnetic field. A detector of internally- 
reflected Cherenkov radiation (DIRC) is used for charged 
particle identification. Kaons are identified with a neural 
network based on the likelihood ratios calculated from 
dE/dx measurements in the SVT and DCH, and from 
the observed pattern of Cherenkov light in the DIRC. 
A finely-segmented CsI(Tl) electromagnetic calorimeter 
(EMC) is used to detect photons and neutral hadrons, 
and to identify electrons. Electron candidates are re- 
quired to have a ratio of EMC energy to track momen- 
tum, an EMC cluster shape, DCH dE/Ax, and DIRC 
Cherenkov angle all consistent with the electron hypoth- 
esis. The instrumented flux return (IFR) contains re- 
sistive plate chambers for muon and long-lived neutral 
hadron identification. Muon candidates are required to 
have IFR hits located along the extrapolated DCH track. 



an IFR penetration length, and an energy deposit in the 
EMC consistent with the muon hypothesis. 



III. DATA SAMPLES 

The data used in this analysis were recorded with the 
BABAR detector at the PEP-II storage ring in the period 
October 1999 to December 2000. The total integrated 
luminosity of the data set is 20.6 fb~^ collected at the 
T(4S') resonance and 2.6 fb~^ collected about 40 MeV 
below the T(4S') (off-resonance data). The corresponding 
number of produced BB pairs is 23 million. 

Samples of Monte-Carlo simulated BB and cc events, 
generated with a GEANT3 ||ll| detector simulation, are 
analyzed through the same analysis chain as the data to 
check for biases in the extracted physics parameters and 
are also used to develop models for describing physics 
and detector resolution effects. However, the values of 
the parameters used in these models are determined with 
data. The equivalent luminosity of this simulated sample 
is approximately equal to that of the data for BB events 
and about half that of data for cc events. In addition, 
we generate signal Monte Carlo samples in which one 
neutral B meson in every event decays to D*~ t^vi, with 
D*^ D^ir^ , and the other neutral B meson decays to 
any final state (l2j. The then decays to one of the 
four final states reconstructed in this analysis (described 
in the next section). The equivalent luminosity of the 
simulated signal samples is between 2 and 8 times that 
of the data, depending on the D° decay mode. 



IV. EVENT SELECTION AND 
CHARACTERIZATION 

We select events containing a fully-reconstructed D*^ 
and an identified oppositely-charged electron or muon. 
This D*~£'^ pair is then required to pass kinematic cuts 
that enhance the contribution of semileptonic B'^ — > 
D*^t^vt decays. In addition to the signal sample, we 
select several control samples that are used to character- 
ize the main sources of background. 

We define the following classification of the sources of 
signal and background that we expect to contribute to 
this sample. The nomenclature shown in italics will be 
used throughout this paper to define signal and all pos- 
sible types of background. Events are classified accord- 
ing to the D*~ candidate reconstruction status and the 
source of the lepton candidate. 

1. Events with a correctly reconstructed D*~ candi- 
date: 

(a) Events that originate from BB events: 

i. Events with a correctly identified lepton 
candidate: 
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A. Signal - B° ^ D*-£+vt (X) decays, 
where the D*~ and lepton originate 
from a common point. {X) indi- 
cates the possibihty of one or more 
pions or photons from the direct de- 
cay of the parent B or from the de- 
cay of short-hved intermediate reso- 
nances (radially- and orbitally-excited 
D states). 

B. Uncorrelated-lepton background - 
events in which the lepton does 
not come from the primary B 
decay that produced the D*^: 
{B D*-X, other B -> 1+ X) or 
{B-^ D*-X, X ^e+Y). 

C. Charged B background - B^ — s- 
D*-i+i^eX. 

ii. Fake-lepton background - events with a 
misidentified lepton candidate. 

(b) Continuum background - cc D*~ X. 

2. Combinatorial- D* feac^rownrf - events with a mis- 
reconstructed D*^ candidate. 



A. Lepton candidates 

Lepton candidates are defined as tracks with momen- 
tum greater than 1.2 GeV/c in the T(45) rest frame. For 
the D*~e+ samples, the electron candidate passes selec- 
tion criteria with a corresponding electron identification 
efficiency of about 90% and hadron misidentification less 
than 0.2%. For the D*~ samples, the muon candidate 
passes selection criteria with a corresponding muon iden- 
tification efficiency of about 70% and hadron misidenti- 
fication between 2% and 3%. The particle identification 
criteria in BABAR are described in detail elsewhere [ p^ . 
A sample enriched in fake-lepton background is also se- 
lected, where D*~£^ candidates are accepted if the lep- 
ton fails both electron and muon selection criteria looser 
than those required for lepton candidates. This sample 
is used to determine the fraction and At distribution of 
the fake-lepton background. 

B. candidates 

D*~ candidates are selected in the decay mode D*^ 
D^Tr~ . The 7?° candidate is reconstructed in the modes 
K~^TT~ , K~^Tr~TT^TT~ , K^Tr~TT^ and K^tt^tt~ The daugh- 
ters of the decay are selected according to the follow- 
ing definitions, tt*^ candidates are reconstructed from two 
photons with energy greater than 30 MeV each, and an 
invariant mass between 119.2 and 150.0 MeV/c^ and a 
total energy greater then 200 MeV. The mass of the 
photon pair is constrained to the 7r° mass and the pho- 
ton pair is kept as a 7r° candidate if the probabil- 
ity of the fit is greater than 1%. candidates are 



reconstructed from a pair of charged particles with in- 
variant mass within 15 MeV/c^ of the mass. The 
pair of tracks is retained as a K'^ candidate if the 
probability that the two tracks form a common vertex is 
greater than 1%. Charged kaon candidates satisfy loose 
kaon criteria fl^ for the K'^n^ mode and tighter cri- 
teria for the if+TT^TT+TT^ and X+tt^tt" modes. For the 
iir+7r~7r° and K'^tt^tt^ modes, a likelihood is calculated 
as the square of the decay amplitude in the Dalitz plot 
for the three-body candidate, based on measured am- 
plitudes and phases |l4j. The candidate is retained if 
the likelihood is greater than 10% of its maximum value 
across the Dalitz plot. This criterion rejects about 95% 
(97%) of uniform background and has a signal efficiency 
of about 62% (48%) for the K+tt-tt° (K'^tt+tt-) mode if 
the real signal is described by the results in Ref. JT^ . 

I?" candidates in the X+vr^, if +7r^7r+7r^, and 
ifgTT+TT^ modes (if+7r~7r° mode) are selected if they 
have an invariant mass within 17 MeV/c^ (34 MeV/c^) 
of the -D'^ mass. The invariant mass of the daughters 
is constrained to the mass and the tracks are con- 
strained to a common vertex in a simultaneous fit. The 
13° candidate is retained if the probability of the fit is 
greater than 0.1%. 

The low-momentum pion candidates for the D*^ 
D'^Tr~ decay are selected with total momentum less than 
450 MeV/c in the T(4S') rest frame and momentum trans- 
verse to the beamline greater than 50 MeV/c. The mo- 
mentum of the D*~ candidate in the T(4S') rest frame is 
required to be between 0.5 and 2.5 GeV/c. The require- 
ments on the momenta of the low-momentum pion and 
D*^ candidates retain essentially all signal events and 
reject higher momentum D*~ from continuum events. 
D*^ candidates are rejected if | cos 6'^|j^jjgj | > 0.85, where 
^thrust is the angle between the thrust axis of the I3*~£+ 
candidate and the thrust axis of the remaining charged 
and neutral particles in the event. The distribution of 
I cos 6**^^^.^,^,^ I is peaked at 1 for jet-like continuum events 
and is fiat for more spherical BB events. 

D*~ candidates are retained if m{D*) — m{D^) is 
less than 165 MeV/c^, where m{D*) is the candidate 
D'-^tt^ mass calculated with the candidate D° mass con- 
strained to the true D'^ mass, m(Z3°). Note that the 
m{D*) — m{D^) distribution has a kinematic threshold 
at the mass of the tt", and a peak at 145.5 MeV/c'^ with 
a resolution of 1 MeV/c^ or better. We have retained the 
sideband of the m{D*) — m{D'^) distribution for studies 
of combinatorial- 13* background. 



C. D*-e+ candidates 

D*^ l'^ candidates are retained if the following criteria 
are met: the x^ probability of the fit of the lepton, 7r~, 
and I?" candidates to a common vertex is greater than 
1%; the decay point of Stag is determined from at least 
two tracks; the fit that determines the distance Az be- 
tween the two B decays along the beamline converges; 
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the time between decays (At) calculated from is less 
than 18 ps; and the calculated error on At (iTAt) is less 
than 1.8 ps. See Sec. ^ for details on the determination 
of the decay point of Stag and the calculation of At and 

CTAt- 

We define two angular quantities for each D*^£^ 
candidate to classify them into a sample enriched in 
Vi signal events, and a sample enriched in 
uncorrelated-lepton background events. The first angle is 
6d-* i, the angle between the D*~ and lepton candidate in 
the T(4S') rest frame. The second is 9b,D'1, the inferred 
angle between the direction of the and the vector sum 
of the D*~ and lepton candidate momenta, calculated in 
the T(45') rest frame. Since we do not know the direc- 
tion of the S", we calculate the cosine of 6b,D'1 from the 
following equation, in which we assume that the only B 
decay particle missed in the reconstruction is a massless 
neutrino: 



cos^ 



D't 



2|pb||pz 



(4) 



All quantities in Eq. ^ arc defined in the T(45') rest 
frame. The energy and the magnitude of the momen- 
tum of the B are calculated from the e+e~ center-of- 
mass energy and the B^ mass. For true B^ D*^ t^Vf^ 
events, cosOb^d-i lies in the physical region [— 1,+1], 
except for detector resolution effects. Backgrounds lie 
inside and outside the range [— 1,-fl]. We also calcu- 
late the same angle with the lepton momentum direc- 
tion reflected through the origin in the T(4S') rest frame: 
9b,d*{-i)- This angle is used to select samples enriched 
in uncorrelated-lepton background. 

A sample enhanced in B^ D*~t^vg signal events 
(called the opposite-side sample) is composed of 
candidates with cosOo-'i < and \ cos9B,D'e\ < 1-1. 
Samples are defined for lepton candidates that satisfy the 
criteria for an electron, a muon and a fake-lepton. The 
first two samples are the signal samples, and the latter is 
the fake-lepton control sample. 

An additional background control sample, representa- 
tive of the uncorrelated-lepton background and called the 
same-side sample, is composed of D*~i'^ candidates sat- 
isfying cosOo'e > and | cos6'b.d*(-^)| < l-l- We use 
cos9b^d-'{-i) rather than cosObm'I because, in Monte 
Carlo simulation, the distribution oicos6B^D*(-i) in this 
control sample is similar to the distribution of cos 9b, D*i 
for uncorrelated-lepton background in the signal sample, 
whereas the distribution of cos 9b, D'i in the background 
control sample is systematically different. 



D. Signal and background subsamples 

Approximately 68, 000 candidates pass the above se- 
lection criteria. These candidates are distributed over 
two signal samples and ten background control samples 
defined by the following characteristics: 



1. whether the data were recorded on or off the T(4S') 
resonance (two choices); 

2. whether the candidate lepton is same-side or 
opposite-side to the D*~ candidate (two choices); 

3. whether the lepton candidate passes the criteria 
for an electron, a muon, or a fake lepton (three 
choices). 

The signal samples are the electron and muon samples in 
the opposite-side, on-resonance data. 

The combinatorial- £** background can be distin- 
guished from events with a real D*^ in a plot of the 
mass difference m{D*) — m{D°). The m{D*) — ■m{D^) 
distributions for the samples of signal events (opposite- 
side _D*^e+ and candidates in on-resonance data) 
are shown as data points in Fig. for (a) electron can- 
didates and (b) muon candidates. The contributions 
of the three types of background that contain a real 
D*^ (continuum, uncorrelated-lepton, and fake-lepton, 
together called the peaking background), except for the 
charged B background, are also shown in the plots. The 
m{D*) — m{D'^) distributions for five background con- 
trol samples used for determining the background levels 
in the signal sample are shown as data points in Fig. |^: 
opposite-side (a) D*~e'^ and (b) D*~ candidates in 
off- resonance data; same-side (c) D*~e+ and (d) D*" fi"^ 
candidates in on-resonance data; (e) opposite-side D*~- 
fake-lepton candidates in on-resonance data. The re- 
maining five background control samples are useful for 
determining the background levels in the first five con- 
trol samples. 

Each of the 12 samples described above are further 
divided into 30 subsamples according to the following 
characteristics that affect the m{D*) — m{D'~') or At dis- 
tributions. 

1. The vr^ from the D*^ decay reconstructed in the 
SVT only, or in the SVT and DCH (two choices): 
The ra{D*) — m{D'^) resolution is worse when the 
7r~ is reconstructed only in the SVT. 

2. The D'^ candidate reconstructed in the mode 
K^TT~ or if+7r~7r° or {K^tt~tt^tt~ or K^tt^-k') 
(three choices): The level of contamination from 
combinatorial-!?* background and the m(D*) — 
m{D^) resolution depend on the D° decay mode. 

3. The 6-tagging information used for the other B (five 
choices; see Sec. Vl): The level of contamination 
from each type of background and the At resolution 
parameters depend on the tagging information. 

This allows subdivision into 360 samples. In the un- 
binned maximum likelihood fits to the m^D*) - m{D°) 
and (At, a At) distributions, individual fit parameters 
are shared among different sets of subsamples based on 
physics motivation and observations from the data. 

We fit the m{D*) — m{D'^) distributions to determine 
signal and background probabilities for each of the 360 
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FIG. 1: m{D*) — m{D'^) distribution for events passing all 
selection criteria for — > D*~£'^ue, with (a) an electron or 
(b) a muon candidate. The points correspond to the data. 
The curve is the result of a simultaneous unbinned maximum 
likelihood fit to this sample of events and a number of back- 
ground control samples. The shaded distributions correspond 
to the four types of background (BG) described in the text. 
The charged B background is not shown separately. 



subsamples. The peak due to real D*^ candidates is 
modeled by the sum of two Gaussian distributions; the 
mean and variance of both the Gaussian distributions, as 
well as the relative normalization of the two Gaussians, 
are free parameters in the fit. Wc model the shape of the 
combinatorial-D* background with the function 



1 

N 



1 — cxp 



6m — 



Cl 



(5) 



where Sm = m{D*) — m{D'^), N is, a. normalization con- 
stant, m^- is the mass of the tt^, and ci and C2 are 
free parameters in the fit. An initial fit is performed to 
determine the shape parameters describing the peak and 
combinatorial- £)* background. Separate values of the five 
parameters describing the shape of the peak are used for 
the six subsamples defined by (1) whether the tt" candi- 
date is tracked in the SVT only or in both the SVT and 
DCH (two choices), and (2) the three types of decay 
modes. Each of these six groups that use separate peak 
parameters is further subdivided into twelve subgroups 
that each uses a different set of the two combinatorial-D* 
shape parameters but the same set of peak parameters. 
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FIG. 2: m{D*) — m{D^) distribution for events passing all 
selection criteria in background control samples; opposite- 
side (a) D*~ and (b) D*~ fi'^ candidates in off-resonance 
data; same-side (c) D*~ and (d) D*~ fi^ candidates in on- 
resonance data; (e) opposite-side D*~-fake-lepton candidates 
in on-resonance data. The points correspond to the data. The 
curve is the result of a simultaneous unbinned maximum like- 
lihood fit to this sample of events, the signal sample, and a 
number of other background control samples. The shaded 
distributions correspond to the four types of background de- 
scribed in the text. The charged B background is not shown 
separately. 



Ten of these twelve subgroups are defined by the five 
tagging categories for the large signal samples and for 
the fake-lepton control samples, in on-resonance data. 
The other two subgroups are defined as same-side, on- 
resonance samples and all off-resonance samples. 

Once the peak and combinatorial-D* shape parame- 
ters have been determined, we fix the shape parameters 
and determine the peak and combinatorial- 13* yields in 
each of the 360 subsamples with an unbinned extended 
maximum-likelihood fit. 

The total peak yields in the signal sample and each 
background control sample are then used to determine 
the amount of true signal and each type of peaking back- 
ground in the m{D*) — m{D°) peak of each sample as 
follows. 

1. Continuum background - For each subsample in 
on-resonance data, the peak yield of the corre- 
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spending subsample in off-resonance data is scaled 
by the relative integrated luminosity for on- and 
off-resonance data, to determine the continuum- 
background yields in on-resonance data. 

2. Fake-lepton background - Particle identification 
efficiencies and misidentification probabilities for 
the electron, muon, and fake-lepton selection cri- 
teria are measured in separate data samples as a 
function of laboratory momentum, polar angle, and 
azimuthal angle, for true electrons, muons, pions, 
kaons, and protons. B^B^ and B^B~ Monte Carlo 
simulations are used to determine the measured 
laboratory momentum, polar angle, and azimuthal 
angle distributions for true electrons, muons, pions, 
kaons and protons that pass all selection criteria 
for D*^£^ candidates, except the lepton (or fake- 
lepton) identification criteria. These distributions 
are combined with the measured particle identifica- 
tion efficiencies and misidentification probabilities 
to determine the momentum- and angle-weighted 
probabilities for a true lepton or true hadron to 
pass the criteria for a lepton or a fake lepton in each 
of the D*~£^ signal and background control sam- 
ples. We then use these efficiencies and misidenti- 
fication probabilities, and the observed number of 
electron, muon and fake-lepton candidates in each 
subsample in data, after removing the continuum 
background contribution, to determine the number 
of true leptons and fake leptons (hadrons) in each 
control sample. 

3. Uncorrelated-lepton background - The relative ef- 
ficiencies for signal and uncorrelated-lepton events 
to pass the criteria for same-side and opposite-side 
samples are calculated from Monte Carlo simula- 
tion. These efficiencies and the m{D*) — m{D^) 
peak yields, after removing the continuum and 
fake-lepton background contributions, are used 
to determine the number of uncorrelated-lepton 
events in each subsample. 

The peak yields and continuum, fake-lepton, and 
uncorrelated-lepton fractions of the peak yield, as well as 
the combinatorial- _D* fraction of all events in a m{D*) — 
m{D^) signal window, arc shown in Table || for the signal 
and background control samples in on-resonance data. 
The peak yields include the peaking backgrounds. The 
signal window is defined as (143-148) MeV/c^ for the cal- 
culation of the combinatorial- _D* background fractions. 
Table || shows the peak yields and the combinatorial-!?* 
background fractions for different divisions of the signal 
sample (opposite-side lepton candidates in on-resonance 
data) . This table demonstrates that the background lev- 
els vary significantly among subgroups of the signal sam- 
ple. 

Finally, we use the calculated fractions and fitted 
shapes of the background sources in each control sample 
to estimate the probability of each candidate to be sig- 



nal or each type of background (combinatorial-!)*, con- 
tinuum, fake-lepton, or uncorrelated-lepton) when we fit 
the (At, CTAt) distribution to determine the lifetime and 
mixing parameters. We take advantage of the fact that 
charged and neutral B decays have different decay-time 
distributions (because the charged B does not mix) to 
determine the fraction of charged B background events 
in the fit to (At, aAt)- 



V. DECAY-TIME MEASUREMENT 

The decay-time difference At between B decays is de- 
termined from the measured separation Az = — ztag 
along the z axis between the D*^t^ vertex position 
{zD'i) and the flavor-tagging decay -Btag vertex position 
(ztag)- This measured Az is converted into At with the 
use of the T(4S') boost, determined from the beam ener- 
gies for each run. Since we cannot reconstruct the direc- 
tion of the B meson for each event, we use the approx- 
imation At « Az/{Pjc). Without detector resolution 
effects, this approximation has a bias that depends on 
the sum of the proper decay times (ti +t2) of the two B 
mesons and their direction in the T(4S') rest frame [ p^ . 
Neither of these quantities can be measured because the 
T(4S') production point is not known and the momen- 
tum of the B is not fully reconstructed due to a missing 
neutrino. After integrating over ti -t- 12 and the B meson 
direction, the mean and RMS of the bias are and 0.2 ps, 
respectively. 

The momentum and position vectors of the D", 7r~, 
and lepton candidates, and the average position of the 
e'^e~ interaction point (called the beam spot) in the 
plane transverse to the beam are used in a constrained 
fit to determine the position of the D*~i~^ vertex. The 
beam-spot constraint is about 100 /xm in the horizontal 
direction and 30 /xm in the vertical direction, correspond- 
ing to the RMS size of the beam in the horizontal direc- 
tion and the approximate transverse flight path of the 
B in the vertical direction. The beam-spot constraint 
improves the resolution on zo'i by about 20% in Monte 
Carlo simulation; the RMS spread on the difference be- 
tween the measured and true position of the D*^£+ ver- 
tex is about 70 /xm (0.4 ps). 

We determine the position of the Stag vertex from all 
tracks in the event except the daughters of the D*^£+ 
candidate, using K'^ —> tt+tt" and A — > pTr~ candidates 
in place of their daughter tracks, and excluding tracks 
that are consistent with photon conversions. The same 
beam-spot constraint applied to the -B_d*£ vertex is also 
applied to the -Btag vertex. To reduce the influence of 
charm decay products, which bias the determination of 
the vertex position, tracks with a large contribution to 
the of the vertex fit are iteratively removed until no 
track has a contribution greater than 6 or only one 
track remains. The RMS spread on the difference be- 
tween the measured and true position of the -Btag ver- 
tex in Monte Carlo simulation is about 160 fim (1.0 ps). 
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TABLE I: Peak yields and the fraction of them that are due to continuum, fake-lepton, and uncorrelated-lepton events. Also 
shown is the combinatorial-D* fraction of total events in a m{D*) — m{D'^) signal window for the signal and background 
control samples in on-resonance data. Peak yields include the peaking backgrounds. The signal window for combinatorial-D* 
background fractions is defined as (143-148) MeV/c^. e, fi, and fake indicate the type of lepton candidate: electron, muon or 
fake-lepton. 



Category Peak Yield /cont(%) /fake(%) /unco.(%) ./•comb(%) 

opposite-side 

e 7008 ±91 1.5 ±0.4 0.168 ±0.004 3.1 ±0.4 17.9 ± 0.2 

At 6569 ±88 2.3 ± 0.6 2.67 ±0.07 2.9 ± 0.5 18.4 ± 0.3 

fake 8770 ±108 12.8 ± 1.3 72.4 ±1.8 0.7 ±1.6 31.4 ±0.2 



e 306 ±21 < 5.9* 0.53 ± 0.04 56.9 ± 7.0 34.0 ± 1.3 
H 299 ± 20 5.1 ±3.6 8.9 ± 0.6 48.9 ± 8.0 34.4 ±1.3 
fake 1350 ±45 20.4 ±4.1 74.4 ± 5.4 3.6 ± 7.8 42.6 ± 0.6 

*90% C.L. 



TABLE II: Peak yields and the combinatorial-D* background 
fraction of total events in a m{D*) — 171(0") signal window for 
different divisions of the signal sample (opposite-side lepton 
candidates in on-resonance data). In the first block, the sig- 
nal sample is divided according to the reconstruction status 
of the 7r~ from the D*~ decay; the second block by the D° 
decay mode; and the third block by the 6-tagging informa- 
tion (see Sec. ^). The signal window for combinatorial-D* 
background fractions is defined as (143-148) MeV/c^. 



Category 



SVT only 
DCH & SVT 

SVT only 
DCH & SVT 



Ktt 

KtYTTTV & KgTVTV 

Ktvk" 



Peak Yield 



5427 ± 81 
1581 ± 41 

5053 ± 78 
1517 ±41 



2623 ± 53 
2219 ± 54 
2166 ±51 



/'comb (%) 



19.5 ±0.3 
11.8 ±0.4 

20.3 ±0.3 
11.1 ±0.4 



7.0 ±0.3 
28.6 ±0.5 
16.9 ±0.5 



Therefore, the At resolution is dominated by the z reso- 
lution of the tag vertex position. 

For each event, we calculate the uncertainty on Az 
(fAz) due to uncertainties on the track parameters from 
the SVT and DCH hit resolution and multiple scattering, 
our knowledge of the beam-spot size, and the average B 
flight length in the vertical direction. The calculated un- 
certainty does not account for errors in pattern recogni- 
tion in tracking, errors in associating tracks with the B 
vertices, the eff'ects of misalignment within and between 
the tracking devices, or the error on the approximation 
we use to calculate At from Az. The calculated uncer- 
tainties will also be incorrect if our assumptions for the 
amount of material in the tracking detectors or the beam- 
spot size or position are inaccurate. We use parameters 
in the At resolution model, measured with data, to ac- 
count for uncertainties and biases introduced by these 
effects. 



VI. FLAVOR TAGGING 



Kn 

Kivmv ( 



lepton 

kaon 

NTl 

NT2 

NT3 

lepton 

kaon 

NTl 

NT2 

NT3 



2491 ± 52 
1939 ± 51 
2139 ± 50 



783 ± 29 
2565 ± 55 
630 ± 27 
921 ± 33 
2108 ±51 

746 ± 28 
2393 ± 53 
545 ± 25 
958 ± 34 
1928 ± 49 



7.4 ±0.3 
30.9 ±0.5 
16.1 ±0.4 



8.2 ±0.6 
17.9 ±0.4 

14.3 ±0.8 
20.9 ±0.7 

20.7 ±0.5 

8.3 ±0.6 
18.6 ±0.4 
15.1 ±0.8 

19.4 ±0.7 

21.8 ±0.5 



All tracks in the event, except the daughter tracks of 
the candidate, are used to determine whether the 

Btag decayed as a B° or a This is called flavor tag- 
ging. We use five different types of flavor tag, or tagging 
categories, in this analysis. The first two tagging cate- 
gories rely on the presence of a prompt lepton, or one or 
more charged kaons, in the event. The other three cat- 
egories exploit a variety of inputs with a neural-network 
algorithm. The tagging algorithms are described briefly 
in this section; see Ref. for more details. 

Events are assigned a lepton tag if they contain an 
identified lepton with momentum in the T(45) rest frame 
greater than 1.0 or 1.1 GeV/c for electrons and muons, re- 
spectively, thereby selecting mostly primary leptons from 
the decay of the b quark. If the sum of charges of all iden- 
tified kaons is nonzero, the event is assigned a kaon tag. 
The final three tagging categories are based on the output 
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of a neural network that uses as inputs the momentum 
and charge of the track with the maximum center-of-mass 
momentum, the number of tracks with significant impact 
parameters with respect to the interaction point, and the 
outputs of three other neural networks, trained to iden- 
tify primary leptons, kaons, and low momentum pions 
from D* decays. Depending on the output of the main 
neural network, events are assigned to an NTl (most cer- 
tain), NT2, or NTS (least certain) tagging category. About 
30% of events are in the NTS category, which has a mistag 
rate close to 50%. Therefore, these events are not sensi- 
tive to the mixing frequency, but they increase the sen- 
sitivity to the i?" lifetime. 

Tagging categories are mutually exclusive due to the 
hierarchical use of the tags. Events with a lepton tag 
and no conflicting kaon tag are assigned to the lepton 
category. If no lepton tag exists, but the event has a 
kaon tag, it is assigned to the kaon category. Otherwise 
events are assigned to corresponding neural network cat- 
egories. The mistag rates are free parameters in the final 
fit. The final results are shown in Table [II in Sec. PQ. 



onance, opposite- or same-side lepton, electron or muon 
or fake lepton - indexed by s) that are further subdi- 
vided into 30 subsamples (tagging category, decay 
mode, with or without DCH hits for the pion from the 
D* decay - indexed by c). We maximize the likelihood 



12 30 N{s,c) 

s=l c=l fe=l 



(6) 



where k indexes the N{s, c) events Xk in each of the 360 
subsamples. The probability Ps,c{Smk, Xk ]ff) of observ- 
ing an event (5mk,Xk), where Xk = (At, (TAt,g), is calcu- 
lated as a function of the parameters 



^ //•comb -'comb -"peak ->comb 



f Pkg r pkg J- pkg 
Js,c,l' s,c,2' J s,c,3' 



i-B+ J ys,c,ii ys,c,2' ys,c,3i yc tIc j 



pkg 



^cli\ 



VII. FIT METHOD 



(7) 



We perform an unbinned fit simultaneously to events 
in each of the 12 signal and control samples (on or off res- as 



P,^ciSmk,Xk;rf) = /,';r'^.F™'"''(,5m;p77'^)g^°'"'^(xfc ;g7,r'') + (1 - /,<;r'').FP'=^^<5m spP^^"^) 



r 



(8) 



where dm is the mass difference ra{D*) — m{D'^) de- 
fined earlier. The symbol "comb" in the first term signi- 
fies combinatorial- £)* background. In the second term, 
the symbol "pkg" denotes peaking background and j in- 
dexes the three sources of peaking background (contin- 
uum, fake-lepton and uncorrelated- lepton). In the last 
term, the parameter /^+ describes the charged- i? frac- 
tion in the sample after all other types of background 
are subtracted, and "sig" and "ch" label functions and 
parameters for the signal and charged-i3 background, 
respectively. The charged-B fraction is assumed to be 
identical for all categories. The index g is -1-1 (—1) for 
unmixed (mixed) events. By allowing different effective 
mistag rates for apparently mixed or unmixed events in 
the background functions g'^°™^ and t/P'^^^ ^fJQ accomo- 
date the different levels of backgrounds observed in mixed 
and unmixed samples. Functions labeled with J- describe 
the probability of observing a particular value of 5m while 
functions labeled with Q give probabilities for values of 
At and UAt in category g. Parameters labeled with / 



describe the relative contributions of different types of 
events. Parameters labeled with p describe the shape of 
a 5m distribution, and those labeled with q describe a 
(At, (TAt) shape. The parameters labeled with p and / 
have been determined by a set of fits to m[D*) — m{D^) 
distributions described in Sec. and are kept fixed in 
the fit to (At, (TAt)- 

Note that we make explicit assumptions that the 5m 
peak shape, parameterized by pP"^"^, and the signal and 
charged- i? background (At, cta*) shapes, parameterized 
by (fc'® and q^^ ^ depend only on the subsample index 
c and not on the control sample index s. The first of 
these assumptions is supported by data, and simplifies 
the analysis of peaking background contributions. The 
second assumption reflects our expectation that the At 
distribution of signal and charged- i? background events 
does not depend on whether they are selected in the sig- 
nal sample or appear as a background in a control sample. 

The ultimate aim of the fit is to obtain the lifetime 
and mixing frequency, which by construction are com- 
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mon to all sets of signal parameters q^^^. Most of the 
statistical power for determining these parameters comes 
from the signal sample, although the fake and uncorre- 
lated background control samples also contribute due to 
their signal content (see Table |) . 

We bootstrap the full fit with a sequence of initial fits 
using reduced likelihood functions restricted to a partial 
set of samples, to determine the appropriate parameter- 
ization of the signal resolution function and the back- 
ground At models, and to determine starting values for 
each parameter in the full fit. 

1. We first find a model that describes the At distribu- 
tion for each type of event: signal, combinatorial- 
D* background, and the three types of backgrounds 
that peak in the m{D*) — m{D^) distribution. To 
establish a model, we use Monte Carlo samples that 
have been selected to correspond to only one type 
of signal or background event based on the true 
Monte Carlo information. These samples are used 
to determine the At model and the categories of 
events {e.g., tagging category, fake or real lepton) 
that can share each of the parameters in the model. 
Any subset of parameters can be shared among any 
subset of the 360 subsamples. We choose param- 
eterizations and sharing of parameters that mini- 
mizes the number of different parameters while still 
providing an adequate description of the At distri- 
butions. 

2. We then find the starting values for the background 
parameters by fitting to each of the background- 
enhanced control samples in data, using the model 
(and sharing of parameters) determined in the pre- 
vious step. Since these background control sam- 
ples are not pure, we start with the purest con- 
trol sample (combinatorial-!?* background events 
from the m{D*) — m{D'^) sideband) and move on 
to less pure control samples, always using the mod- 
els established from earlier steps to describe the At 
distribution of the contamination from other back- 
grounds. 

The result of the above two steps is a At model for each 
type of event and a set of starting values for all param- 
eters in the fit. When we do the final fit, we fit all sig- 
nal and control samples simultaneously (approximately 
68,000 events), leaving all parameters in the Q functions 
free in the fit, except for a few parameters that either are 
highly correlated with other parameters or reach their 
physical limits. The total number of parameters that are 
free in the fit is 72. The physics parameters Tgo and Am^ 
were kept hidden until all analysis details and the sys- 
tematic errors were finalized, to minimize experimenter's 
bias. However, statistical errors on the parameters and 
changes in the physics parameters due to changes in the 
analysis were not hidden. 



VIII. SIGNAL At MODEL 

For signal events in a given tagging category c, the 
probability density function (PDF) for At consists of a 
model of the intrinsic time dependence convolved with a 
At resolution function: 

^-s(At,aAt,<?;gT') = 



1 



g-|Att.uc|/r^o [1 + _ 2uj^) C0s(AmrfAtt,.uc)] 



4tbo 

(g)7^((5A^, f7At;C'^), 



(9) 



where 7?, is a resolution function, which can be different 
for different event categories, g is +1 (—1) for unmixed 
(mixed) events, SAt is the residual At — Attruc, and uJc 
is the mistag fraction for category c. To account for an 
observed correlation between the mistag rate and UAt in 
the kaon category (described in Sec. VIII A), we allow 
the mistag rate in the kaon category to vary as a linear 
function of aAt- 



C^kaonf At 



offset 



(10) 



and the offset to 



and allow both the slope 
be free parameters. In addition, we allow the mistag 
fractions for 5° tags and tags to be different. We 
define Aw = w^o —ujgo and uj = (w^o +uj^o)/2, so that 



^BO/BO = ^ ^ ■ 



(11) 



The model for the intrinsic time dependence has 13 pa- 
rameters: LOc and AlOc for each of the five tagging cate- 
gories, ckkaon, Arrid and tbo. 

For the At resolution model, we use the sum of a single 
Gaussian distribution and the same Gaussian convolved 
with a one-sided exponential to describe the core part of 
the resolution function, plus a single Gaussian distribu- 
tion to describe the contribution of "outliers" — events 
in which the reconstruction error SAt is not described by 
the calculated uncertainty aAt- 

7eGExp+G(<5At, aAt; s, K, /, 6°"*, s™*, /™*) 
= fG{SAt;0,.saAt) 

+ (1 - / - /°"')G(u - SAt; 0, saAt) ® Eiu; naAt) 

-h/™*G((5At;6™*,s™*) , 

(12) 

where u is an integration variable in the convolution G (8) 
E. The functions G and E are defined by 

G{x; xo,a) = exp [-{x - xo f/{2a f] (13) 

yzTia 

and 



E{x; a) 



— j a '^■^P i^l^) if 2; < 0, 



ff a; > 0. 



(14) 
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The exponential component is used to accommodate a 
bias due to tracks from charm decays on the Stag side. 

Since the outUer contribution is not expected to be de- 
scribed by the calculated error on each event, the last 
Gaussian term in Eq. ^ does not depend on cta* ■ How- 
ever, in the terms that describe the core of the resolution 
function (the first two terms on the right-hand side of 
Eq. p^ ), the Gaussian width s and the constant k in the 
exponential are scaled by (XAt . The scale factor s is intro- 
duced to accommodate an overall underestimate (s > 1) 
or overestimate (s < 1) of the errors for all events. The 
constant n is introduced to account for residual charm de- 
cay products included in the i3tag vertex; k is scaled by 
CTAt to account for a correlation observed in Monte Carlo 
simulation between the mean of the SAt distribution and 
the measurement error cta*- 

The correlation between SAt and ciAt is due to the 
fact that, in B decays, the vertex error ellipse for the D 
decay products is oriented with its major axis along the 
D flight direction, leading to a correlation between the 
D flight direction and the calculated uncertainty on the 
vertex position in z for the Stag candidate. In addition, 
the flight length of the D in the z direction is correlated 
with its flight direction. Therefore, the bias in the mea- 
sured Stag position due to including D decay products 
is correlated with the D flight direction. Taking into ac- 
count these two correlations, we conclude that D mesons 
that have a flight direction perpendicular to the z axis in 
the laboratory frame will have the best z resolution and 
will introduce the least bias in a measurement of the z 
position of the Stag vertex, while D mesons that travel 
forward in the laboratory will have poorer z resolution 
and will introduce a larger bias in the measurement of 
the Stag vertex. 

The mean and RMS spread of At residual distributions 
in Monte Carlo simulation vary significantly among tag- 
ging categories. We find that we can account for these 
differences by allowing the fraction of core Gaussian, /, 
to be different for each tagging category. In addition, 
we find that the correlations among the three parame- 
ters describing the outlier Gaussian (6°"*, s°"*, /°"*) are 
large and that the outlier parameters are highly corre- 
lated with other resolution parameters. Therefore, we 
fix the outlier bias 5°"* and width s°"', and vary them 
over a wide range to evaluate the systematic uncertainty 
on the physics parameters due to fixing these parameters 
(see Sec. KII). The signal resolution model then has eight 
free parameters: s, k, /°"*, and five fractions /c (one for 
each tagging category c). 

As a cross-check, we use a resolution function that is 
the sum of a narrow and a wide Gaussian distribution. 
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and a third Gaussian to describe outliers: 

7^G+G+G(<5A^, aAub, s, /, b^, s'", 5™*, /™*) 

= fG{SAt;baAt,scrAt) 

+ (!-/- r')G{SAt; b^aAus'^^At) 

+ /°"*G((5At;5™*,s™*) . 
This resolution function has two more parameters than 
7?.GExp+G- It accommodates a bias due to tracks from 
charm decays on the Stag side by allowing the means of 
the Gaussian distributions to be nonzero. 



A. Vertex-tagging correlations 

A correlation dwc/daAt ~ 0.12 ps~-^ is observed be- 
tween the mistag rate and the At resolution for kaon 
tags. This effect is modeled in the resolution function 
for signal as a linear dependence of the mistag rate on 
a At, as shown in Eq. ^ In this section, we describe the 
source of this correlation. 

We find that both the mistag rate for kaon tags and the 
calculated error on At depend inversely on -^/Spf , where 
Pt is the transverse momentum with respect to the z axis 
of tracks from the Stag decay. Correcting for this depen- 
dence of the mistag rate removes most of the correlation 
between the mistag rate and a At- The mistag rate de- 
pendence originates from the kinematics of the physics 
sources for wrong-charge kaons. The three major sources 
of mistagged events in the kaon category are wrong-sign 
mesons from S decays to double charm (b ccs), 
wrong-sign kaons from Z)+ decays, and kaons produced 
directly in S decays. All these sources produce a spec- 
trum of tracks that have smaller \/Spf than S decays 
that produce a correct tag. The CAt dependence origi- 
nates from the 1 /p^ dependence of cr^ for the individual 
contributing tracks. 



IX. At MODELS FOR BACKGROUNDS 

Although the true At and resolution on At are not 
well-defined for background events, we still describe the 
total At model as a "physics model" convolved with a 
"resolution function" . 

The background At physics models we use in this anal- 
ysis are all a linear combination of one or more of the 
following terms, corresponding to prompt, exponential 
decay, and oscillatory distributions: 
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^ph"*(^*true,5) 
^'phys(^*true,5) 
a°|%(Aitrue,5) 



-5(Att™e) [l + 5(l-a;P™*)] , 

J^exp(-|Att™e|/r'^«) [l + 5(l-a;'"-^)] , 

^ exp(-|Ait,ue|/T'^^) [1 + 5(1 - 0;°^^) COs{Am''^Attrne)] , 



(16) 

(17) 
(18) 



r 



where 6 (At) is a (5- function, = +1 for unmixed and —1 
for mixed events, and t^^ and Am^^ are the effective hfe- 
time and mixing frequency for the particular background. 

For backgrounds, we use a resolution function that is 
the sum of a narrow and a wide Gaussian distribution: 

nG+G{SAt,aAub,s,f,b^,s^) 

= fG{SAt; baAt, sa^t) + (1 - f)G{6At; b^a^t, s"<TAt) ■ 

(19) 



A. Combinatorial-f* background 

Events in which the D*~ candidate corresponds to a 

random combination of tracks (called combinatorial-D* 
background) constitute the largest background in the sig- 
nal sample. We use two sets of events to determine 
the appropriate parameterization of the At model for 
combinatorial-!?* background: events in data that are in 
the upper m(D*) ^ m{D'^) sideband (above the peak due 
to real D*~ decays); and events in Monte Carlo simula- 
tion that are identified as combinatorial-D* background, 
based on the true information for the event, in both the 
m{D*) — m{D^) sideband and peak region. The data 
and Monte Carlo At distributions are described well by 
a prompt plus oscillatory term convolved with a double- 
Gaussian resolution function: 



■'Comb 



j^osc^osc /A+ _comb 



...comb ^jOsc 



1 



+ (1 - r^)a^hys(Attrue,5;C^P'-'"*)] <^TIg+G ■ 

(20) 



Approximately 60% of combinatorial-D* background 

events are from B'^B^ events according to Monte Carlo 
simulation. Although the D*~ is not correctly recon- 
structed, the identified lepton is very likely to be a pri- 
mary lepton. The tagging algorithm can still identify the 
flavor of Stag with a reasonable mistag probability, espe- 
cially for the lepton category, and for the kaon category 
if the tracks swapped between the D*~i~^ candidate and 
Stag are pions. Therefore, the combinatorial-!?* back- 
ground also exhibits oscillatory behavior. 

The parameters w?™*, Am'=°'"'^, t"""^^, f, 6"", and s'" 
are shared among all subsamples. The parameters 0;°'*'^, 
/°'''^, b, and s are allowed to be different depending on 
criteria such as tagging category, whether the data were 



recorded on- or off-resonance, whether the candidate lep- 
ton passes real- or fake-lepton criteria, and whether the 
event passes the criteria for same-side or opposite-side 
D*~ and i. The total number of free parameters in the 
combinatorial- D* background At model is 24. 

The relative fraction of B^B'^ and B^B^ events in 
the combinatorial-!?* background depends slightly on 
m{D*) — m{D^). However, no significant dependence of 
the parameters of the At model on m{D*) — m{D°) is 
observed in data or Monte Carlo simulation. The sam- 
ple of events in the m{D*) — m{D^) sideband is used to 
determine the starting values for the parameters in the 
final full fit to all data samples. 

To reduce the total number of free parameters in the 
fit, parameters that describe the shape of the wide Gaus- 
sian (bias and width) are shared between combinatorial- 
D* background and the three types of peaking back- 
ground: continuum, fake-lepton, and uncorrelated- 
lepton. The wide Gaussian fraction is allowed to be dif- 
ferent for each type of background. 



B. Continuum peaking background 

All cc events that have a correctly rccoustriicted D*^ 
are defined as continuum peaking background, indepen- 
dent of whether the associated lepton candidate is a real 
lepton or a fake lepton. The cc Monte Carlo sample and 
off-resonance data are used to identify the appropriate 
At model and sharing of parameters among subsamples. 
The combinatorial-!)* background At model and param- 
eters described in the previous section are used to model 
the combinatorial-!?* background contribution in the off- 
resonance At distribution in data. 

The decay vertex of a real !?*^ from continuum cc pro- 
duction always coincides with the primary vertex. If the 
lepton candidate also originates from the primary ver- 
tex, we can use a prompt physics model convolved with 
a resolution function that can accommodate a bias due to 
tracks from charm decays other than the !?*^ candidate. 
If the lepton candidate is from a charm decay, the mea- 
sured vertices of the D*~£~^ candidate and the remaining 
tracks are both likely to be between the primary vertex 
and the charm vertex; hence the measured Az is likely 
to be very small. Both types of events can be modeled 
with a prompt model convolved with a double-Gaussian 
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resolution function: 



D. Uncorrelated-lepton peeiking background 



gcont ^ gprrnt^^^^^^^^ ^. ^prmt) ^ ^^^^ _ (2I) 



Dependence on the flavor tagging information is included 
to accomodate any differences in the amount of back- 
ground events classified as mixed and unmixed. 

By fitting to the data and Monte Carlo control sam- 
ples with different sharing of parameters across subsets of 
the data, we find that the apparent "mistag fraction" for 
events in the kaon category is significantly different from 
the mistag fraction for other tagging categories. We also 
find that the core Gaussian bias is significantly differ- 
ent for opposite-side and same-side events. We introduce 
separate parameters to accommodate these effects. 

The total muiibcr of parameters used to describe the 
At distribution of continuum peaking background is six. 
The off-resonance control samples in data are used to 
determine starting values for the final full fit to all data 
samples. 



C. Fake-lepton peaking background 

To determine the At model and sharing of parameters 
for the fake-lepton peaking backgrounds, we use B^B^ 
and B~^B~ Monte Carlo events in which the D*^ is cor- 
rectly reconstructed but the lepton candidate is misiden- 
tified. In addition, we use the fake-lepton control sample 
in data. Th(^ combinatorial-D* and continuum peaking 
background At models and parameters described in the 
previous two sections are used to model their contribu- 
tion to the fake-lepton At distribution in data. For this 
study, the contribution of signal is described by the sig- 
nal parameters found for signal events in the Monte Carlo 
simulation. 

Since the fake-lepton peaking background is due to B 
decays in which the fake lepton and the D*^ candidate 
can originate from the same B or different B mesons, 
and the charge of the fake lepton can carry correct flavor 
information of the reconstructed B candidate, we include 
both prompt and oscillatory terms in the At model: 



Qfake 



^osc/^osc /"A J. ^. ^fake A^take , ,osc\ 
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(22) 



We find that the apparent mistag rates for both the 
prompt and mixing terms, and the bias of the core Gaus- 
sian of the resolution function, are different between some 
tagging categories. The total number of parameters used 
to describe the fake-lepton background is 14. The fake- 
lepton control samples in data are used to determine 
starting values for the final full fit to all data samples. 



To determine the At model and sharing of parameters 
for the uncorrelated-lepton peaking backgrounds, we use 
B°B° and B+B' Monte Carlo events in which the £>*" is 
correctly reconstructed but the lepton candidate is from 
the other B in the event or from a secondary decay of 
the same B. In addition, we use the same-side control 
sample in data, which is only about 30% uncorrelated- 
lepton backgroimd in the m(_D*)— m(D*') peak region due 
to significant contributions from combinatorial-!?* back- 
ground and signal. The combinatorial-!?* and peaking 
background At models and parameters described in the 
previous two sections are used to model their contribu- 
tion to the same-side At distribution in data. For this 
initial fit, the contribution of signal is described by the 
signal parameters found for signal events in the Monte 
Carlo simulation. 

Physics and vertex reconstruction considerations sug- 
gest several features of the At distribution for the 
uncorrelated-lepton sample. First, we expect the recon- 
structed At to be systematically smaller than the true 
At value since using a lepton and a D*~ from different 
B decays will generally reduce the separation between 
the reconstructed B^'e and i?tag vertices. We also ex- 
pect that events with small true At will have a higher 
probability of being misreconstructed as an uncorrelated 
lepton candidate because it is more likely that the fit of 
the D*~ and ^ to a common vertex will converge for these 
events. Finally, we expect truly mixed events to have a 
higher fraction of uncorrelated-lepton events because in 
mixed events the charge of the D* is opposite that of 
primary leptons on the tagging side. These expectations 
are confirmed in the Monte Carlo simulation. 

We do not expect the uncorrelated-lepton background 
to exhibit any mixing behavior and none is observed in 
the data or Monte Carlo control samples. We describe 
the At distribution with the sum of a lifetime term and 
a prompt term, convolved with a double-Gaussian reso- 
lution function: 



j'lifez-'life / A j. ^. 
/ yphys(^*true,5';T 

phys 



+ (1 - /"■°)a::hvs(Attrue,5;'^P™*)] ® ^G+G • 

(23) 



The effective mistag rates a;P''™* and w^'^** accommodate 
different fractions of uncorrelated-lepton backgrounds in 
events classified as mixed and unmixed. We find that 
the apparent mistag rate for the lifetime term is different 
between some tagging categories. All other parameters 
are consistent among the different subsamples. The total 
number of parameters used to describe the uncorrelated- 
lepton background is six. The uncorrelated-lepton con- 
trol samples in data are used to determine starting values 
for the final full fit to all data samples. 
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E. Charged B peaking background 

The charged-i3 peaking background is due to decays 
of the type D*£i'(X. Since charged i?'s do not 

exhibit mixing behavior, we use the At and tagging in- 
formation to discriminate charged- i? peaking background 
events from neutral- i? signal events, in the simultaneous 
fit to all samples. We use the same resolution model and 
parameters as for the neutral-B signal since the At res- 
olution is dominated by the ztag resolution and the B 
decay dynamics are very similar. The charged B back- 
ground contribution is described by 



1 



^-=^e-|A-^-l/^.Ml + .9(l-2c.|..)] ^^^^ 

(E)TZ{SAt,aAt;qc) , 

where is the mistag fraction for charged B mesons 
for tagging category c. 

Given that the ratio of the charged B to neutral B life- 
time is close to 1 and the fraction of charged B mesons 
in the peaking sample is small, we do not have sufficient 
sensitivity to distinguish the lifetimes in the fit. We pa- 
rameterize the physics model for the B~^ in terms of the 
lifetime ratio /t^o , and fix this ratio to the Review 
of Particle Properties 2002 world average of 1.083[^. We 
vary the ratio by the error on the world average (±0.017) 
to estimate the correspon ding systematic uncertainties 
on Tgo and Anid (see Sec. XII ). 

In each tagging category, the fit is sensitive to only 
two parameters among the neutral B mistag frac- 

tion {ojb") and the charged B fraction {fB+)- Therefore 
we fix the ratio of mistag rates, ujb+/lobo, to the value 
of the ratio measured with fully reconstructed charged 
and neutral hadronic B decays in data, for each tagging 
category. 



X. FIT RESULTS 

The total number of free parameters in the final fit is 
72: 21 in the signal model, one for the charged B fraction, 
24 in the combinatorial- £)* background model, and 26 in 
peaking background models. The fitted signal At model 



parameters are shown in Table III. 

The statistical correlation coefficient between r^o and 
Arud is p{Amd, Tgo) = —0.22. The global correlation co- 
efficients (the largest correlation between a variable and 
every possible linear combination of other variables) for 
Tgo and Amd, and some of the correlation coefficients 
between r^o or Ato^ and other parameters, are shown in 



Table IV 



Figure || shows the At distributions for unmixed and 
mixed events in a sample in which the probability of each 
event being a signal is higher than a threshold chosen so 
that the sample is 80% pure in signal events. The points 
correspond to data. The curves correspond to the sum of 



TABLE III: Results of full fit to data — signal model and 
resol ution function parameters. A correction, described in 
Sec. XI A, has been applied to Tgo and Amd- The uncertain- 



ties are statistical only. 



parameter 


value 


parameter 


value 


Amd (,ps j 






— U.iiz ± U.Uzo 


TgO (ps) 


1 1-90 +0.024 

— U.Uzd 




-0.023 ±0.019 


/s+ 


0.082 ± 0.029 


S 


1.201 ± 0.063 


^lepton 


0.071 ±0.015 




0.86 ±0.17 


offset 
^kaon 


0.002 ± 0.024 


/"lepton 


0.72 ±0.10 


Qkaon (ps"''") 


0.229 ± 0.036 


./kaon 


0.609 ± 0.088 




0.212 ±0.020 


/nti 


0.69 ±0.13 




0.384 ±0.018 


/nT2 


0.70 ±0.10 




0.456 ±0.012 


/nT3 


0.723 ± 0.078 


AcJiepton 


-0.001 ±0.022 


yout 


0.0027 ± 0.0017 


AWkaon 


-0.024 ±0.015 


(ps) 


-5.000 (fixed) 


Al^UTl 


-0.098 ±0.032 


(ps) 


6.000 (fixed) 



TABLE IV: Global correlation coefiicients for Amd and Tgo 
from the full fit to data and other correlation coefficients for 
pairs of key parameters in the fit. 



Amd global correlation 
TgO global correlation 

p(Amd,Tso) 
p{Amd,fB+) 
p(Tso.Ssig) 



0.74 
0.69 

-0.22 
0.58 
-0.49 
-0.26 



the projections of the appropriate relative amounts of sig- 
nal and background At models for this 80%-pure signal 
sample. Figure ^ shows the time-dependent asymmetry 



A{At) = 



ed 

(At) - Ar,nixcd(At) 



N^nmi^ediAt) + N.^i^^diAt) ' 



(25) 



The unit amplitude for the cosine dependence of A is 
diluted by the mistag probabilities, the experimental At 
resolution, and backgrounds. 

Figure ^ shows the At distributions for unmixed and 
mixed events, and the asymmetry A(At) for data sam- 
ples in which events are selected based on the background 
probabilities such that the sample contains 99.5%-pure 
combinatorial background events (left plots), or 60%- 
pure fake- lepton background events (right plots). The 
observed oscillatory behaviors are expected as explained 
in Sec. 



IX 



Since many parameters in the model are free, it is in- 
teresting to see how the errors on Tgo and Amd, and their 
correlation, change when different parameters are free in 
the fit, or fixed to their best value from the full fit. We 
perform a series of fits, fixing all parameters at the val- 
ues obtained from the default fit, except (a) Anid and 
Tgo, (b) Arud, Tgo, and all mistag fractions in the signal 
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FIG. 3: The At distribution on linear (a, b) and logarithmic (c, d) scale for (a, c) unmixed and (b, d) mixed events in an 
80%-pure signal sample and the projection of the fit results. Each event in this sample has a probability of being a signal higher 
than a threshold chosen so that the sample is 80% pure in signal events. The shaded area shows the background contribution 
to the distributions. 
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FIG. 4: The asymmetry plot for mixed and unmixed events 
in an 80%-pure signal sample and the projection of the fit re- 
sults. Errors on the data points are computed by considering 
the binomial probabilities for observing different numbers of 
mixed and unmixed events while preserving the total number. 



background parameters adds a very small contribution to 
the error. The contribution from the charged B fraction 
and mistag fractions to the r^o error is negligible. On 
the other hand, the charged B fraction changes the error 
on Arud the most. The contributions from floating the 
mistag fractions, resolution functions, and background 
At models are relatively small. 

We also check the statistical errors on data by measur- 
ing the increase in negative log likelihood in data in the 
two-dimensional (r^o , Arn^) space in the vicinity of the 
minimum of the negative log likelihood. We find that the 
positive error on Tgo is about 6% larger than that deter- 
mined by the fitting program, whereas the other errors 
are the same as those determined by the fit. To take this 
into account, we increase the positive statistical error on 
tbo by 6%. 



XI. VALIDATION AND CROSS CHECKS 



model, (c) Am^, r^o, and (d) Am^, t^o, and 
all mistag fractions in the signal model, (e) all parame- 
ters in the signal At model. The one-sigma error ellipses 
for these fits and for the default fit are shown in Fig. |^. 

We can see that the error on r^o changes very little 
until we float the signal resolution function. Floating the 



In Sec. XI A , we describe several tests of the fitting pro- 
cedure that were performed with both fast parameterized 
Mon te Ca rlo simulations and full detector simulations. In 
Sec. XI B| , we give the results of performing cross-checks 
on data, including fitting to different subsamples of the 
data and fitting with variations to the standard fit. 
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FIG. 5: The At distributions for (a, d) unmixed and (b, e) mixed events, and (c, f) the asymmetry plot in a 99.5%-pure 
combinatorial-D* sample (a, b and c) and in a 60%-pure _D*~-fake-lepton event sample (d, e and f). Events are selected based 
on the background probabilities, such that the sample contains 99.5%-pure combinatorial-/)* events, or 60%-pure D*~-iake- 
lepton background events. The projection of the fit results is overlayed on top of the data points. Errors on the data points 
in the asymmetry plots are computed by considering the binomial probabilities for observing different numbers of mixed and 
unmixed events while preserving the total number. 



Tests of fitting procedure witli Monte Carlo 
simulations 



A test of the fitting procedure is performed with fast 
parameterized Monte Carlo simulations, where 87 exper- 
iments are generated with signal and background control 
sample sizes and compositions corresponding to that ob- 
tained from the full likelihood fit to data. The mistag 
rates and At distributions are generated according to the 
model used in the likelihood fit. The full fit is then per- 
formed on each of these experiments. We find no statis- 
tically significant bias in the average values of r^o and 
Arud for the 87 fits. The RMS spread in the distribution 
of results is consistent with the mean statistical error 
from the fits and the statistical error on the results in 
data, for both tbo and Am^. We find that 20% of the 
fits result in a value of the negative log likelihood that is 
smaller (better) than that found in data. 

We also fit two types of Monte Carlo samples that in- 
clude full detector simulation: pure signal and signal plus 
background. To check whether the selection criteria in- 
troduce any bias in the lifetime or mixing frequency, we 



fit the signal physics model to the true lifetime distribu- 
tion, using true tagging information, for a large sample 
of signal Monte Carlo events that pass all selection cri- 
teria. We also fit the measured At distribution, using 
measured tagging informat ion, w ith the complete signal 
At model described in Sec. VIII . We find no statistically 
significant bias in the values of Tgo or Anid extracted in 
these fits. 

The B°B°, B+B-, and cc Monte Carlo samples that 
provide simulated background events along with signal 
events are much smaller than the pure signal Monte Carlo 
samples. In addition, they are not much larger than the 
data samples. In order to increase the statistical sensitiv- 
ity to any bias introduced when the background samples 
are added to the fit, we compare the values of Tgo and 
Amd from the fit to signal plus background events, and 
pure signal events from the same sample. We find that 
when background is added, the value of r^n increases by 
(0.022 ± 0.009) ps and the value of Anid increases by 
(0.020 ± 0.005) ps^^, where the error is the difference 
in quadrature between the statistical errors from the fit 
with and without background. We correct our final re- 
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-0.02 -0.01 0.01 0.02 

5Amd(ps ) 



TABLE V: Summary of systematic uncertainties on the two 
physics parameters, Tgo and ArUd- 



Source 


5(Amd) (ps"-") 


S{TBa) (ps) 


Selection and fit bias 


0.0123 


0.0178 


z scale 


0.0020 


0.0060 


boost 


0.0005 


0.0015 


Beam spot position 


0.0010 


0.0050 


SVT alignment 


0.0030 


0.0056 


Background / signal prob. 


0.0029 


0.0032 


Background At models 


0.0012 


0.0063 


Fixed B+/B° lifetime ratio 


0.0003 


0.0019 


Fixed / mistag ratio 


0.0001 


0.0003 


Fixed signal outlier shape 


0.0010 


0.0054 


Signal resolution model 


0.0009 


0.0034 


Total systematic error 


0.013 


0.022 



FIG. 6: Comparison of one-sigma error ellipses in the Am^- 
Tgo plane for fits in which different sets of parameters are 
free. From the innermost to the outermost ellipse, the float- 
ing parameters are (Am^, t^o), (Arritj, r^o, mistag fractions), 
(Amd, r^o, (Amd, Tgo, fg+, mistag fractions), all sig- 

nal At parameters, and the default fit (72 floating parame- 
ters). 



suits in data for these biases, which are each roughly 
the same size as the statistical error on the results in 
data. We conservatively apply a systematic uncertainty 
on this bias equal to the full statistical error on the mea- 
sured result in Monte Carlo simulation with background: 
±0.018 ps for tbo and ±0.012 ps^^ for Am^. 



B. 



Cross-checks in data 



We perform the full maximum-likelihood fit on differ- 
ent subsets of the data and find no statistically significant 
difference in the results for different subsets. The fit is 
performed on datasets divided according to tagging cat- 
egory, &-quark flavor of the D*~t^vi candidate, 
&-quark flavor of the tagging B, and decay mode. We 
also vary the range of At over which we perform the fit 
(maximum value of \A.t\ equal to 10, 14, and 18 ps), and 
decrease the maximum allowed value of cta* from 1.8 ps 
to 1.4 ps. Again, we do not find statistically significant 
changes in the values of r^o or Am^. 



XII. SYSTEMATIC UNCERTAINTIES 

We estimate systematic uncertainties on the parame- 
ters and Amji with studies performed on both data 
and Monte Carlo samples, and obtain the results sum- 
marized in Table 



The largest source of systematic uncertainty on both 
parameters is the limited statistical precision for deter- 
mining the bias due to the fit procedure (in particular, 
the background modeling) with Monte Carlo events. We 
assign the statistical errors of a full fit to Monte Carlo 
samples including bac kgrou nd to estimate this systematic 
uncertainty. See Sec. XI A for more details. 

The calculation of the decay-time difference At for each 
event assumes a nominal detector z-scale, PEP-II boost, 
vertical beam-spot position, and SVT internal alignment. 
The PEP-II boost has an uncertainty of 0.1% based on 
our knowledge of the beam energies. The z-scale uncer- 
tainty is determined by reconstructing protons scattered 
from the beam pipe and comparing the measured beam 
pipe dimensions with the optical survey data. The z- 
scale uncertainty is less than 0.4%. We shift the vertical 
beam-spot position by up to 80 /im, or vary the posi- 
tion randomly with a Gaussian distribution with a width 
of up to 80 /xm, and assign the variation in the fitted 
parameters as a systematic uncertainty. The systematic 
uncertainty due to residual errors in SVT internal align- 
ment is estimated by reprocessing the simulated sample 
with different internal alignment errors. We assign the 
shift of the fitted parameters as a systematic uncertainty. 

The modeling of the ra{D*) — m{D^) distribution de- 
termines the probability we assign for each event to be 
due to signal. We estimate the uncertainty due to the sig- 
nal probability calculations by repeating the full fit using 
an ensemble of different signal and background param- 
eters for the m{D*) — m{D^) distributions, varied ran- 
domly according to the measured statistical uncertainties 
and correlations between the parameters. We assign the 
spread in each of the resulting fitted physics parameters 
as the systematic uncertainty. 

The modeling of the background At distribution affects 
the expected background contributions to the sample. 
The systematic uncertainty due to the assumed back- 
ground At distributions is estimated as the shift in the 
fitted parameters when the model for the largest back- 
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ground (due to combinatorial-D* events) is replaced by 
the sum of a prompt term and a lifetime term. 

The model of the charged B background assumes fixed 
B'^ /B^ ratios for the mistag rates and lifetimes. We vary 
the mistag ratio by the uncertainty determined from sep- 
arate fits to hadronic B decays. We vary the lifetime ra- 
tio by the statistical uncertainty on the world average Q . 
The resulting change in the fitted physics parameters is 
assigned as a systematic uncertainty. 

The final category of systematic uncertainties is due to 
assumptions about the resolution model for signal events. 
We largely avoid assumptions by fioating many param- 
eters to describe the resolution simultaneously with the 
parameters of interest. However, two sources of system- 
atic uncertainty remain: the shape of the outlier con- 
tribution, which cannot be determined from data alone, 
and the assumed parameterization of the resolution for 
non-outlier events. We study the sensitivity to the outlier 
shape by repeating the full fit with outlier Gaussian func- 
tions of different means and widths. The mean is varied 
between — 1 ps and —10 ps, and the width is varied from 
4 ps to 12 ps. We assign the spread of the resulting fitted 
values as a systematic uncertainty. We estimate the un- 
certainty due to the assumed resolution parameterization 
by repeating the full fit with a triple-Gaussian resolution 
model and assigning the shift in the fitted values as the 
uncertainty. 

The total systematic uncertainty on tbo is 0.022 ps and 
on Ato^ is 0.013 ps~^. 



Both the lifetime and mixing frequency have combined 
statistical and systematic uncertainties that are compa- 
rable to those of the most precise previously-published 
experimental measurements The results are con- 

sistent with the world average measurements of Tgo = 
(1.542 ± 0.016) ps and Am^ = (0.489 ± 0.008) ps^^ §. 

This analysis demonstrates the feasibility of measuring 
the B^ lifetime and mixing frequency simultaneously at 
B Factory experiments, realizing the advantages of bet- 
ter determinations of the At resolution function and the 
amount of B^ background. All background fractions. 
At resolution parameters for signal and background, and 
mistag fractions are determined from the data. The life- 
time is most correlated with the At resolution parameters 
for signal, while the mixing frequency is most correlated 
with the B^ background fraction. The largest system- 
atic uncertainty on both parameters is the limited sta- 
tistical precision for determining any bias due to the fit 
procedure (in particular, the background modeling) with 
Monte Carlo simulation. 

Both the statistical and systematic uncertainties on 
these parameters can be reduced with the larger data and 
Monte Carlo simulation samples already available at the 
B Factories. Other physics parameters, such as the dif- 
ference in decay rates of the neutral B mass eigenstates, 
can also be included in a simultaneous fit in future data 
samples. 



XIII. SUMMARY 

We use a sample of approximately 14,000 exclusively 
reconstructed B^ D*~t'^vg signal events to simulta- 
neously measure the B'^ lifetime tbo and oscillation fre- 
quency Am^. We also use samples of events enhanced 
in the major types of backgrounds. The lifetime and 
oscillation frequency are determined with an unbinned 
maximum-likelihood fit that uses, for each event, the 
measured difference in decay times of the two B mesons 
and its uncertainty, the signal and background proba- 
bilities, and 6-quark tagging information for the other 
B. In addition to the lifetime and oscillation frequency, 
we extract the parameters describing the signal At res- 
olution function, the background At models, the mistag 
fractions, and the B^ background fraction, in the simul- 
taneous fit to signal and background samples. The results 
for the physics parameters are 

Tbo = (1.523lg;^23 ± 0.022) ps 

and 

Awd = (0.492 ± 0.018 ± 0.013) ps^^ 

The statistical correlation coefficient between r^o and 
Amd is -0.22. 
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